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Abstract 

In this work we investigate charge transport in one-dimensional (ID) insula- 
tors via semi-classical and perturbative renormalization group (RG) methods. 
We consider the problem of electron-electron, electron-phonon and electron- 
two-level system interactions. We show that non-linear collective modes such 
as polarons and solitons are reponsible for transport. We find a new excita- 
tion in the Mott insulator: the polaronic soliton. We discuss the differences 
between band and Mott insulators in terms of their spin spectrum and obtain 
the charge and spin gaps in each one of these systems. We show that electron- 
electron interactions provide strong renormalizations of the energy scales in 
the problem. 
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I. INTRODUCTION 



It is with joy that we celebrate Michael Pollak's birthday. One of us (A. H.C.N.) had 
the honor of inheriting Mike's office after his retirement. In this office we had long and 
exciting conversations about many aspects of condensed matter physics and especially one 
of Mike's early interests: DNA electronics. It turns out that Mike worked for a few years as 
an experimentalist and he was probably one of the first researchers to look into the problem 
of DNA stability with condensed matter experimental tools ffl. These were probably one of 
the first systematic studies of DNA molecules from the condensed matter point of view. In 
fact, Mike went ahead and also worked on some theoretical aspects of DNA ||- Besides being 
a complete scientist, theorist and experimentalist, Mike's interest in DNA was far ahead of 
his own time. Mike's broad vision of science can also be seen in his beautiful and important 
work on the insulator-metal transition [§]. Perhaps because of his knowledge and experience, 
Mike was always convinced that DNA was not a simple band insulator but a more exotic 
Mott insulator. At this point in time it is not clear if a DNA molecule is an ordinary band 
insulator or the more complex Mott insulator. We believe that the evidences that DNA 
is an insulator are quite strong: DNA is transparent to light and transport measurements 
have shown a charge gap of the order of a few volts @,[J. From the theoretical point of 
view the question that remains is: is DNA a band insulator or a Mott insulator? While 
naive electron counting and LDA calculations |J seem to indicate that DNA is a band 
insulator with a very narrow bandwidth, there are some experimental indications that DNA 
is not an ordinary band insulator. The first comes from the fact that antiferromagnetic 
excitations have been observed in DNA J7| (see our discussion of magnetic fluctuations in 
Mott insulators below) and DNA seems to be able to conduct supercurrents when connected 
to superconducting electrodes via a proximity effect ||. These effects cannot be observed 
in ordinary insulators where the gap to charge excitations produces a vanishing density of 
states at the chemical potential. So, maybe, once again Mike was right and DNA is hiding 
a few surprises. However, the field of DNA electronics is still in its infancy. Many more 
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experiments are needed in order to understand the real nature of electron propagation in 
DNA and the theoretical effort is intense P-jE^. It is with an eye to these possibilities that 
in the last few years we have been trying to understand the problem of conduction in ID 
insulators, which is the main theme of this paper. We are thankful to Mike for pointing the 
way. 

Band insulators cannot transport charge at low temperatures because of a charge gap 
in the excitation spectrum. When doped with acceptors and/or donors impurities, a band 
insulator can conduct charge as in the case of semiconductors. But the periodicity of the 
lattice is broken and impurity bands can only be formed at high density of dopants. There is, 
therefore, a competition between the disorder (which localizes the carriers), and the overlap 
between the impurities wave-function (which delocalizes the carriers). In band insulators 
electron-electron interactions are supposed to be weak and are not usually discussed. While 
in three dimensions this is certainly true, ID systems are clearly strongly interacting because 
electrons cannot avoid each other in their motion, in other words, interactions are enhanced 
in ID because of phase space constraints. In this paper we are going to discuss the case of 
clean insulators in one dimension that can only transport charge via non-linear excitations. 
By their nature these non-linear excitations only exist if the electronic system interacts with 
itself or with another set of degrees of freedom that can provide feedback effects and hence 
nonlinearity. 

The classical example of a non-linear excitation in a insulator is a polaron: the dressing 
of a single electron in the conduction band by phonons. Because of non-linear effects the 
polaron problem is equivalent to the self-trapping of the electron and the creation of a bound 
state below the conduction band (very much like a donor state in a semiconductor). The 
polaron problem has a long history and it is a well-known non-linear problem |13| . Another 
famous problem is the propagation of solitons in polymers like polyacetylene. While poly- 
acetylene is an insulator because of the Peierls mechanism that leads to the doubling of the 
unit cell, doping can produce lattice-soliton states in the middle of the gap that in principle 
can carry charge [f[4] , |l5|| . As in the case of band insulators, the Peierls mechanism does not 
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require any electron-electron interactions and usually the lattice distortion is assumed to 
be static so polaronic band effects are not usually discussed. In systems where insulating 
behavior is driven by electron-electron interactions, that is in the case of the so-called Mott 



insulators, charge-solitons are known to exist ||T6| , |iT|| . Usually charge-solitons are again 
only discussed in the context of a static lattice without any polaronic effects. These are 
collective excitations associated with the electronic charge density and, like polarons and 
lattice-solitons, they are energetically costly since these are topological defects of the field 
theory that describes them in the first place. As we will show there are many similarities 
between polaron and soliton conduction and in general these two types of excitations play 
an important role in ID systems. 

In this paper we generalize the problem of soliton formation in ID insulators by taking 
into account not only the electron-electron interaction (which is unavoidable in these sys- 
tems) but also the electron-phonon interaction and the interaction between electrons and 
two-level systems. Our interest in discussing the interaction between two-level systems and 
electrons arises from three main sources. The first is the electron interaction with hydrogen 
bonds in biological systems such as DNA [T8j , the second is the interaction between elec- 



trons and dangling bonds of atoms in ID nano-structures pi, and finally two level systems 



appear in problems of magnetic moments in ID metals - the so-called Kondo chain p0| . In 
all cases the electronic scattering by localized two-level systems can lead to new effects such 
as unitary scattering and Kondo resonance. Thus, two-level systems behave like pseudo- 
spins. As is expected, the level of complexity of the systems we are going to discuss is quite 
high but, as we are going to show, the interplay between the electrons with these degrees 
of freedom leads to new effects that cannot be observed in their absence. As we are going 
to see the energy for creation of non-linear excitations is greatly reduced by feedback effects 



21] and this might be at the core of the problem of electron propagation in DNA molecules 
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The paper is organized as follows: in Section |IJ we describe the model; in Section [TTJ we 
describe the problem of a ID non-interacting band insulator (this is a somewhat idealized 



situation but will illustrate very well the main concepts in polaron physics and how they 
apply in one-dimension); in Section [IV] we discuss the problem of non-linear transport in a 
Mott insulator and compare with its non- interacting counterpart; in Section |VJ we show that 
along the so-called Luther-Emery line, where the interacting electronic problem is exactly 
solvable, that the electron-phonon coupling reduces to the problem of polaron dressing of 
the soliton excitations (the polaronic soliton); Section [VI] contains our conclusions and a 
discussion of the interplay between electron-electron interaction and the lattice potential. 



II. DESCRIPTION OF THE MODEL 

Our system of interest can be broken into three different pieces: electrons, phonons and 
pseudo-spins (two-level systems). The electrons are described in term of creation, ip^x), and 
annihilation, tp a (x), operators at position x with spin a (o =|, j.) and obey anti-commutation 
relations: j^o- ( x ) j VC' ( x ')} — S (x — x') S^y. In the absence of interactions the Hamiltonian 
describing the electron motion is simply: 

Ho = J dxj: {|^^ + nxwiwM*)} , (2-i) 

where m is the electron mass and V(x) is the lattice potential. If V(x) is periodic then a 
gap opens up in the spectrum and if the number of electrons per unit cell is 2 the problem 
is described band insulator. 

While the above description is satisfactory for a system where the electronic wave- 
function is extended (and therefore the system is well-described by an electron gas), in 
systems where the coupling with the ions is strong a tight-binding description is usually a 
better starting point. In fact this kind of description is the starting point for the description 
of the Mott insulator. The electron Hamiltonian including the electron-electron interactions 
can be generically written as: 

H e = -t e C+ a C l+1<a + n iAni,i + V ^ n it(7 n i+1 y , (2.2) 

i,<7 i i^o^u' 
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where Cj )(7 {C^ a ) is the annihilation (creation) electron operator at unit cell % with spin 
a and = Cf a C ij(T is the electron number operator. t e is the hopping energy of the 
electron between different sites and U and V describe the on-site and the nearest neighbor 
electron-electron interactions, respectively. 

In the non- interacting limit the electron energy is Ej~ = —2t e cos(ka), where a is the 
lattice spacing. For a finite density of electrons all the states up to the Fermi energy, Ep, 
are filled generating two Fermi points defined by E±k F = Ep, where kp = nn/2 is the Fermi 
momentum (n is the number of electrons per unit of length). Close to the Fermi points 
the spectrum is Ek F + q — Ep = vpq where vp = 2t e a sin (kpa), is the Fermi velocity. The 
linearization of the spectrum close to the Fermi points can also be translated into operator 
language. We rewrite the electron operator as the product of a rapidly varying part ^e l±kpX ) 
and a slowly varying part: 

Ms) = i>RAxV kFX + ^{xY~ lkFX > (2-3) 

where ^r j(T (x) (ipL,a( x )) creates a right (left) moving electron in the system. Because of the 
chiral nature of these excitations they can be described in terms of bosonic operators : 

i> {RtL) M = -^e±*^(W).-W . (2.4) 

In turn, the bosonic modes (j)(nL),o can be described in terms of new fields, (p a and 9 a , as 
4>{R,L),a{ x ) — 0o-( a; ) T 9o-( x )- The bosonic fields can then be rewritten in terms of charge and 
spin bosonic modes: 

$P,<T = 7! (0T ± 0i) 

e Pj(T = ± 6 {) . (2.5) 

The procedure of writing electron operators in terms of bosons is called bosonization. It can 
be shown that the charge density operator p (x) is related to $ p (a?) by 

V2 d% (x) 

p(x = . 2.6 

7T OX 
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The non-interacting Hamiltonian and all the scattering process in the ID system can now be 
rewritten in terms of bosons. There are three main types of scattering in ID: forward scat- 
tering (with small momentum transfer), backscattering (with 2k f momentum transfer) and 
Umklapp scattering (with 4k f momentum transfer). While the forward and backscattering 
do not require the lattice to participate, the Umklapp scattering requires, by momentum 
conservation, that 4kF = G where G = 2n/a is the reciprocal lattice vector. Because of this 
constraint the Umklapp scattering is only important at commensurate filling factors and is 
responsible for the opening of the Mott gap in the half- filled system (n = 1/a). 

The model ( p.2|) can be bosonized using the fields <3> p (x) and $ CT (x). The Hamiltonian 
of the problem breaks into two main parts, H = H s + H c , describing spin H s and charge H c 
23|: 



H s = H a + 
H c = H + 



20i 



(2vra) 

2<?3 



(27TOT 

In the above Hamiltonian we have 



dx cos 
dx cos 



V8§p (x) + 4k F x 



(2.7) 



H v = — dx 
2tt 



U V K V ) 7T 2 n^ (x) + — ^ 



u v ( d$ v (x) 



V = p,(T 



(2.8) 



K v \ dx J 

describing a gapless Luttinger liquid involving only the forward scattering processes. Here, 
u v are the velocity of the bosonic fields, K v are the Luttinger parameters, and (x) is 
the momentum field operator that is canonically conjugate to (x) : [$„ (x) , H v > (x r )] = 
i6 v yS (x — x'). 



In (|2.7| ) the coupling constants g\ and represent the backward scattering and the 
Umklapp process. One can relate the parameters in the bosonized Hamiltonian with the 
ones in Q by g|: 

Va 

UpK p = u a K a = v F 

7T 



Up/ K p = v F + 
u a /K a = v F - 



{U + 5V)a 

71 

(U-V)a 

7T 



01 = Hi = u - 2V. . (2.9) 
a a 

This relation between the coupling constants in the bosonic theory and the interactions in 
the electron problem allow us to change freely from one representation to the other. Notice 
that the bosonization procedure is not a good starting point for the problem of the band 
insulator because n = 2/ a and therefore Vp = 0, implying that we cannot define right and 
left moving particles (this is just the trivial statement that when the band is full there is no 
electronic motion). 

Observe that the bosonic fields in (|2.7| ) separate completely into the spin and charge 
parts. This is called spin-charge separation and it is a generic property of one dimensional 
systems. Namely, while an isolated electron carries both spin and charge, in the ID system 
the electron decays into the bosonic collective modes described above that carry spin and 
charge separately. Since we are only considering the properties of insulators, the spin degrees 



of freedom decouple from the problem and are completely described by H s in ( 2.7|) . From this 
result we can immediately draw some conclusions about the nature of the spin excitations 
in these systems. 

The simplest way to understand the importance of these terms is to perform a pertur- 
bative renormalization group (RG) calculation by assuming that the non-linear terms in 
( 2/7 ) are small compared to the Luttinger liquid terms in (|2.8|) . We can study how these 
non-linear terms affect the physics of the Luttinger liquid. The RG procedure is defined 
by introducing a cut-off A in momentum space so that all the momenta, k, are defined in 
—A < k < A. We now trace all degrees of freedom in the region of high momenta around 
the cut-off, say, A/6 < \k\ < A (b > 1), and rescale the frequencies and momenta so that the 
Luttinger liquid Hamiltonian (|2.8| ) becomes invariant under rescaling [uo — > buo and k — > bk) . 
As we rescale the system to lower energies the coupling constants change with the scale. 
This change is given by the RG equations. For the backscattering problem it can be shown 



that [23 



^M = 2(l-K a ) 9l (e), (2.10) 



where I = ln(l/A). Notice that for K a > 1, the backscattering term scales to zero and 
therefore it is irrelevant under the RG. The spin degrees of freedom are completely described 
by the Luttinger liquid Hamiltonian ( |2.8| ). The spin excitations are gapless and propagate 
with velocity v a . However, when K a < 1, the backscattering term grows under the RG, 
implying that the non-linear terms are becoming stronger at low energies: the backscattering 
is relevant. The relevance of this term can be easily understood if we consider only the non- 
linear part of the Hamiltonian: 

R i = ^r 1 ^ I dx cos{Vs{^{x)) . (2.11) 
l\jxa) J 

If g 1 — > oo at long length scales we see that the coefficient of the cosine grows and becomes 
much bigger than the energy scale of the Luttinger liquid in Q2.8| ). If we now ignore com- 
pletely the Luttinger liquid part of the Hamiltonian (as the RG procedure tells us we should 
do), we see that the energy of the problem is minimized when the field $ CT (a;) = for g± < 



or $> a (x) = y7r/8 when g% > 0. Thus, in the ground state the field is uniform. Let us allow 
the field $ CT to fluctuate around its mimimum value. Expanding Hj up to second order in 
$„- we find 



« + r4 . (2 - 12) 

2(irar (nar 1 



where L is the size of the system. The first term in ( 2.12|) is the ground state energy of the 



uniform state. Adding now the Luttinger liquid part of the Hamiltonian, Eq. ( |2.8| ), we find 
that the action that describes the system is: 



\gi\L r dg r dw_ 



2(na) 2 ' J 2n J 2tt 



$ CT (g,u;)| 2 , (2.13) 



2, V ° 2 , 2 \9l\ 



2K a v a " 2K a ' (ira)\ 
where we have Fourier transformed the fields in the last line. Notice that the excitation 
spectrum has changed from u(q) = ±v a \q\ to 



u(q) = ±v a \Jq 2 + m 2 



m = —JK a \g 1 \v (7 . (2.14) 
7ra v 



This result indicates that a gap of energy 



A s = v rT m = 2\^^ (2.15) 
V ixa 

opens in the spectrum at q = (k = kp in terms of the original electrons). We can define a 
correlation lenght £ a — 1/m — v a /Ao which tells us that the correlations decay exponentially 
over £ CT and that for distances larger than £ CT no correlations are possible. Thus, the relevance 
of the operators under the RG implies that the low lying spin excitations of the system are 
gapped (in contrast to the Luttinger liquid). This is the so-called spin gap and is a sign of 
Cooper pairing in the system since for K a < 1 the interactions become effectively attractive. 
To understand the origin of this result we notice that when two electrons form a Cooper 
pair in a singlet state it costs energy to break the pair and make a spin excitation. 
Let us focus on the charge part of the Hamiltonian that can be written as: 

2" 

+ 

[2-ko] 

The problem described by (|2.16|) is well understood and it can be shown that when 4k pa ^ 2-7T 
the Umklapp term is irrelevant and the problem is described by a gapless Luttinger liquid. 
However, when the system becomes commensurate with the lattice (n = 1/a) the Umklapp 
term in ( |2.16| ) is important. In fact, an RG calculation gives [ p3| : 

dg 3 (£) 



H e = ±- [dx (u p K p ) 7r 2 nJ + ^ f^p] + -^-2 [ dx cos \V8$ p + Ak F x] . (2. 16) 
2n J K p \ ox I (2iraY J 



di 



2(1 - K p )g 3 (i) . (2.17) 



In complete analogy with the backscattering problem, this result shows that a gap opens 
in the charge spectrum (the so-called Mott gap) if K p < 1. The Umklapp term is irrel- 
evant when Kp > 1. Thus, a correlated insulator has two major requirements. Namely, 
the electron-electron interactions have to be repulsive and the charge density has to be 
commensurate with the lattice. 

In fact the calculation of the gap in ( |2.15| ) is incorrect because we have not taken into 
account the fluctuations in the problem by expanding the fields close to its minimum. We 
can learn more about the actual value of the Mott gap by looking at the RG equation (|2.17|) 
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more closely. That RG predicted that for K p < 1 the system scales to strong coupling and 
#3 grows without bound. The maximum value g% can attain in our theory is the largest 
energy scale, that is, the electron bandwidth, W ~ 4i e . Beyond this energy scale the 
perturbative RG calculation breaksdown. In our linearized theory it is easy to see that 
W ps vpA , where A ~ l/o is the bare cut-off in the theory. When gs(£) becomes of the 
order of W the variable I reaches its maximum value £* = ln(Ao/A*). Here A*, which has 
dimensions of inverse length, becomes of the order of the inverse of correlation length, that 
is, A* ps l/£ p ps A/vp. However, accordingly to (|2.17|) we must have 

1 



2(1 - K p , 



Hg(n/gsm 



HW/A) ps l _ R HWj\g z \) , (2.18) 
which can be solved for A as 

Am = W (2.19) 

showing that the gap is proportional to \g^\ 1 ^ 2i ' 1 ~ Kp ^ instead of l^l 1 ^ 2 as predicted by our 
naive expansion (|2.15| ) . In fact the square root behavior of the gap with the coupling constant 
is only obtained when K p = which is equivalent to the strong coupling regime where the 
interaction becomes of the order of the bandwidth of the original problem (see (|2.9f) ). The 
exact value of Am can be obtained by the mapping, via a Jordan- Wigner transformation, 



the Umklapp problem into the anisotropic Heisenberg model |25]. This model can be solved 
exactly by Bethe ansatz and the gap calculated without approximations |16|- It is found 
that the power law dependence of the Mott gap with the coupling constant in ( |2.19| ) is exact. 
This result shows the importance of the RG calculation. Observe that the bandwidth W 
which appears in ( |2.19| ) is a non-universal number that cannot be obtained from the RG. 



By comparing ( |2.19| ) with ( |2.15[ ), however, we find W ps Avp/iira). 



The other important components of the problem are the phonons. These can be of two 
types: acoustic or optical. Since we are only interested in the low energy physics of the ID 



system we consider only the acoustic modes that can be described by the Hamiltonian |26 
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H ph = - I d.v 



2 



n|(x) 2 (d<t> L {x)\ 2 



[2.20) 



Pl 3 \ dx J 

where Pl is the mass density of the system, c s is the sound velocity, and IIl (x) is the 
phonon momentum operator that is canonically conjugated to the phonon field 4>l{x): 
[0 L (x),n L {x')]=iM (x-x'). 

The electron-phonon coupling assumed here is of the deformation potential type, which 
is appropriate for acoustic phonons: 

d<p L (x) 



H e - P = j p J dx p (x) 



dx 



/d(j)[j (xj 
dx V> f (x) ip (x) — — — 

lp V2 [dx d%{x)d4> L { x ) 



tt J dx dx 

Another component in our problem is the presence of pseudo-spins (or two-level systems) 
that couple to the electronic charge degrees of freedom but not to the spins (this is why we 
term them pseudo-spins). In a real system these pseudo-spins can be dangling bonds of atoms 
in a nano-structure or hydrogen bonds in DNA. Because of their nature the pseudospins are 
assumed to be non-interacting and are described by a Hamiltonian: 

H t = l -Jdx[-s(j z (x) +ta x (x)], (2.22) 



where a z and a x are Pauli matrices that obey: [a z (x), a x (x')] = 2ia y (x)5 (x — x'). In (|2.22|) 



e > is the energy difference between the ground state with pseudo-spin j and the first 
excited state with pseudo-spin j. For generality we also allow for the quantum tunneling of 
the pseudo-spin between the two possible configurations and t is the energy associated with 
this tunneling. 

As mentioned previously, the coupling between the pseudo-spins and the charge degrees 
of freedom has essentially the same structure as the electron-phonon coupling: 



H e-t = fjdx^ (X)V(Z) [o*fr)-l] 



dx — [a z {x) - 1] . (2.23) 



\/2tx J dx 
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The rationale for this coupling is the fact that electron polarization around a pseudo-spin 
can lead to its local reorientation. This is the same reasoning behind the electron-phonon 
coupling except that in the case of the pseudo-spin the reorientation process is quantized. 

In order to simplify the problem we project out the pseudo-spin variables in terms of a 
new field 9{x) by writing the generic form for the pseudo-spin wavefunction, 



and rewriting the Hamiltonian as: 



^ cos6(x) ^ 



sm8(x) 



(2.24) 



/ 



H t -> ($\H t \V) = \j dx [-ecos29(x)+tsm29(x)} 
H e _ t -> (9\H e - t \V) = -7* J dx^{x)ij){x) sin 2 (x) . 



(2.25) 



Now, all the variables in the problem are given in terms of fields that are continuous and 
smoothly varying in space. 



III. POLARONS IN BAND INSULATORS 

In this section we review the problem of polarons in band insulators described by the 
Hamiltonian ( |2.1| ). The problem reduces to a valence band which is full and a conduction 
band that is empty. In this discussed previously, the bosonization procedure is not 

a good starting point since the charge density is such that the Fermi velocity vanishes and 
therefore the linearization around the Fermi point is a bad approximation. The problem of 
bosonization and band insulators for weak periodic potentials is discussed in Section [V]]. In 
this section we disregard electron-electron interactions and consider the problem of a single 
electron (hole) at the bottom (top) of the conduction (valence) band. The Hamiltonian 
for the electron in the conduction band is identical to the one in (|2.1|) except that we can 
set V = and replace the electron bare mass m by its effective band mass m c . The full 
Hamiltonian of the problem is written as H p = H + H p h + H t + H e _ p + H e _ t where: 
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H 
H p h 



2m r 



dx 



dijfi dijj 
dx dx 



- I dx 

2 

1 



a (d<j> L 



H t = dx [-£ cos 20 (a;) + t sin 20(e)] 
H e _ t = — 74 / dx^ (x)ip(x) sin 2 (x) . 



The generating function of the problem can be written as: 

Z = J Dij*DipD(j) L D6 e -iSdtjdxCW,i>,<l>L,e] ? 

where ip* and ip are Grassmann variables and the Lagrangean of the problem is: 



r -t,(i* d ^ i 9 ^ 



+ - [e cos 26 - t sin 29] - 



h 2 



2m r , \ dx J \dx 

L 



PL 

2 



Of 



dx 



1p- 



dx 



— 74 sin 6{x) 



(3.1) 



(3.2) 



(3.3) 



Although the above Lagrangean is of relative complexity we are going to study the problem 
in the semiclassical limit by letting h — > 0. In this limit, as one can see directly from ( |3.2| ), 
the problem is dominated by the saddle point equations: 



t cos 20 = -e sin 29 + % a |^| 2 sin 26> 



ih 



Pl- 



d 2 



d^ 
~dt 

L 



h 2 d 2 ip 
2m c d 2 x 



+ [%- 



dt 2 



C 2 sPL 



d 2 <p L d\ 



dx 

i 1 2 



jt sin 2 9 ) ip(x) 



dx 2 



+ 1 P - 



The first equation in (|3.4|) gives: 

tan20(o;) 



dx 



t 



(3.4) 



(3.5) 



which shows that the pseudo-spins can rearrange their orientation due to the presence of 
the electrons. In fact if we define a position dependent pseudo-spin energy: 



six) = e - 74 \mx) 



(3.6) 
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we see that the electrons decrease the effective pseudo-spin energy allowing for the tunneling 
of the pseudo-spins. In other words, the electrons delocalize the pseudo-spins. In what 
follows we are going to consider two regimes. In the first regime we have t 3> e(x) and 
therefore the pseudospins delocalize, in the second regime t -C e the renormalization of the 
pseudo-spin energy is not enough to delocalize the pseudo-spins. As we are going to see, 
in both limits the physics is very similar and therefore the final behavior of the polaron is 
essentially independent of these limits. The expression ( |3.5|) can also be written as: 

sin 2 6(x) = - ; . (3.7) 

2 2 Je*(x)+t* 



For t ^> e we have 

sin 2 6(x)^±-( lt mx)\ 2 -e + t) , (3.8) 



2t 

while for t <C e we have 



2 , t { _ l././_„N |2 , e 



sm 2 6(x)^ — ^ t mx)\ z + -) . (3.9) 



Substituting Q or (Q back into O we find: 

P^ = c sPL^+l P ^ r , (3.10) 

where the values of D t and E t are giving by 

1 2 

E t-\( l -l)lt (3.11) 



2 V t 
for £ ^> £, and 



t£ i7 t (3-12) 
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for i < e. The form of the equations in (|3.10|) is similar to the ones studied by Davydov 
p7| in the problem of propagations of solitonic waves in biological systems by taken into 
account electrons and phonons. This should not be surprising since they describe non- 
interacting electrons interacting with harmonic excitations. We should stress, however, that 
our equations also include the pseudo-spin variables that were not considered in ref. |27| . 

The solution of the problem can now be obtained by noticing that both the electron 
field, ip, and the phonon field, 0l, have the form of a traveling wave. Using the Galilean 
invariance of the problem we write: 



ip (x, t) = 0o (x T vt) exp 
i> L {x, t) = (j) L (x =f vt) = (j) L (A) 



/ 7-1 \ 

— (mvx — E t) 
n 



(3.13) 



where A = x =F vt. Here Eq is the binding energy of the electron and v is the speed of the 
traveling wave. After this simple transformation we find: 



h 2 rf 2 c 
2m r d\ 2 



+ 7 P 



PL [V ~C 



Combining both equations, we have 



d(j) L , 
dX 00 = 


E - 


A d(j) L 




v dX ' 





+ E t 



+ D t <p 3 



(3.14) 



h 2 d 2 (j) 

2m r dX 2 



E Q - — + E t 



+ 



PL (C 2 S - V 2 ) 



*8, 



(3.15) 



which is the non-linear Schrodinger equation supporting the propagation of a soliton in the 
system. This soliton is nothing but a polaron formed of the polarization of the lattice and 
the pseudo-spins. It is simple to show that for one electron (or hole) the solution is: 

g (x - vt) 



00 OM) = -^sech 



(3.16) 



where 



D t + 



PL (C 2 S ~ V 2 ) 



(3.17) 



and the binding energy is 
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En 



2^2 



(3.18) 



By comparing this result with the one found in reference [2l| where pseudo-spins are not 
considered, we notice the binding energy is lowered by E t . Thus, the soliton becomes more 
stable due to the electron coupling with two-level systems. Moreover, since the size of the 
polaron is of order £ « 1/g we see that the coupling to the pseudo-spins also decreases the 
size of the polaron by a factor of D t . The relative lattice displacement is given by ( |3.14j ) 



<j>L(x,t) 



1p 



2Pl (c 2 s 



tanh 



g(x- vt) 



(3.19) 



IV. POLARONIC SOLITON IN MOTT INSULATORS 



In the previous section we studied the problem of polaron formation in a band insulator. 
In this section we will study the problem of non-linear excitations in a Mott insulator where 
the gap is generated by correlations, not by the interaction with the lattice. Since this is 
an interacting problem it is clearly more involved. However, because we are considering a 
system where the electronic density is commensurate with the lattice and the band is not 
full, we can use the tools of bosonization to study it. The methods will be, however, very 
similar to the ones used in the previous section. This will show that the polaronic soliton 
formation in the Mott insulator is very similar to the problem of polaron formation in a 
band insulator. 

The basic Hamiltonian of the problem is given by: H m = H ee + H p h + H t + H e _ p + H e _ t 
where: 

2" 



2tt 



dx 



K n \ dx 



ph 



dx 



+ 



2^3 



(27rct) 



dx cos 



H t = ^Jdx [-e cos 26 (x) + t sin 29{x)\ 



H, 



e—p 



1p / dx 



d§ p (x) d<p L (x) 



dx 



dx 
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r d<$> p (x) . 2 
H e - t = ~1t J dx — ^- — sin 6 [x 



(4.1) 



Once again the problem can be discussed in terms of the generating functional for the 
fields: 



Z = J D<& p D(j) L Dd e %JdtfdxC[$ P ,<p L ,e] ^ 



(4.2) 



where the Lagrangean of the system can then be written as 



u p K p 



dt 



u, 



'd<$> f 



K„ \ dx 



2g 3 



+ i [ecos2 ,_ isiD291 _ 7 ^ 



The equations of motion are: 



+ 7r 



(2tt« 
<9$ 



■ cos 



2 



dx 



dx 



sin 2 9 . 



(4.3) 



t cos 29 = -e sin 29 + 7t^ sin 29 

dx 



d 2 (j) L dsin 2 9 



Up-wKp dt 2 nKp dx 2 



[Tea 



dx 2 



7r 



dx 



Pl- 



d 2 <j> L 2 d 2 <p L d 2 <$> r 



(4.4) 



dt 2 s dx 2 dx 2 

The first equation in ( |4.4|) can be solved exactly as in ( p3|) and the effective local pseudo- 
spin energy is replaced by 

d$ n 



e\x) = s - 7 t - 



dx 



(4.5) 



We consider the same limits as before, namely t ^> e and t < e. The equations of motion 
for the fields are: 



1 d 2 <$> r 



d 2 $ P , V2g 3 . 



UpTiKp dt 2 ixKp dx 2 (txol) 



PL 



d\ 



dx 2 



dx 2 



dt 2 



dx 2 



(4.6) 



This set of equations is very different from the ones obtained in the case of the band insulator 
that is described by the Davydov's equations ( |3.1(J| ). In fact these equations contain non- 
perturbative effects that cannot be described by the the non-interacting electron picture. 
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Using once again the Galilean invariance of the solutions we can write $ p (x, t) = 
$ p (x =F vt) = $ p (A) and 4>l ( x , t) = 4>l ( x T vt) = 4>l (A), where A = x =F vt. The equations 
of motion reduce to: 

+ D t\—r7T = - ^sm(V8$p +7 P - 



UpTxKp nK p 7 dA 2 (naY v 'WA 2 

^ u " C J ^xT = 7 P -nF- (4-7) 



c/A 2 IP dX 2 

Combining both equations and replacing $ p with $ = V8& p , the equation of motion for the 
charge field becomes: 

d 2 $ 



dX 2 

where sign(x) = +1 (—1) if x > (x < 0) is a sign function, 



sz£/n(/x)sin($(A)), (4.8) 



A = ±^=A (4.9) 



with 



(7ra) 



3 



Pl (c 2 - t; 2 ) Wp 71 "-^ 



(4.10) 



Notice that ( |4.8| ) is a sine-Gordon equation that has soliton solutions ||21|| : 
$ p (x, t) = ±V2 tan" 1 fe^ (a:-V * )N ) , for /i > 



= iv^tan" 1 (eVM {x + JL=, for fi < 0. (4.11) 

The electronic density that is related to $ p by equation ( |2.6| ) and is localized in space as 
in the polaron case ( |3.16|) with a width of order £ rs 1 / \/Ta*T- I n the present case, however, 
the soliton properties are heavily dressed by the lattice and pseudo-spin. Therefore we will 
term this excitation a polaronic soliton. 



V. THE POLARONIC SOLITON AT THE LUTHER-EMERY LINE 



In the previous section we described the physics of the polaronic soliton formation in 
terms of the bosonized theory. Since the connection between the electrons and boson fields 
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is not linear, it is not straightforward to interpret the formation of solitons in the bosonic 
language as formation of solitons in the fermionic language. Moreover, it is not clear what 
is the dressing process that changes an electronic soliton into a polaronic soliton. To make 
the connection more explicit we are going to study the problem of electron-phonon coupling 
along the so-called Luther-Emery line. Along this line the bosonic degrees of freedom can 
be re-fermionized. To do so we re-introduce the backscattering terms into the Hamiltonian 
and the parameters in the Hamiltonian are chosen so that [ [2^ , [29| ,j30[l : 

9i ~ 2 92 = 6 

TTVf + #4 5 ' 

where 

9l = 9± = U + 2 V. 
a a 

This is equivalent to choosing K p = 1/2. The charge part of the problem can be exactly 
diagonalized via an unitary transformation and the final Hamiltonian reads : 

Hle = ^E k (p[ pk C lpk - Cl pk C 2p k) , 
k 

where C\ pk {C\ k ) is the transformed annihilation (creation) operator for momentum k for 
left (right) moving electrons. The energy spectrum is given by 



^=4 



(5.1) 



Obviously the spectrum has a gap of magnitude Am = J^-. This is the Mott gap. Notice 
that because K p = 1/2 this value of the gap is in agreement with the RG result in ( |2.19| ). 
We now add electron-phonon interaction ( |2.21| ) into the Hamiltonian as 



H LE . ph = J2E k (cl pk C lpk - C\ pk C 2pk ) - ^ E \/ P P (P) (ap + a ~v) > ( 5 -2) 

where a], (or a p ) is the phonon creation(or annihilation) operator, M is the mass of a unit 
cell, and 

p(p) = E (p\ pk+p Ci pk + Cl pk+p C 2p k) (5.3) 

k 
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is the transformed fermionic density. Notice that the unitary tranformation that diagonalizes 
the problem does not modify the form of the density operator. What this result shows is 
that at the Luther-Emery problem becomes identical to the polaron problem discussed in 
Section [[II]. Since the operators C produce the solitonic excitations in the electronic system 
it becomes clear that the polaronic soliton is formed in two steps: (1) the electron-electron 
interaction generates collective non-linear modes; (2) these non-linear modes are dressed by 
the lattice fluctuations. 

Let us make this description more concrete and assume that a soliton is created in the 
upper band of the Mott insulator. We can expand (|5.1|) for small q as: 

E >^ + &' <«> 

where 

m* = f^V , Am „ (5.5) 

V 4 J 2(ttv f + g 4 ) 2 V ; 

is the effective mass of the carrier. Let us assume that the electron-phonon coupling is weak 

and compute the change in energy of the ground state due to this coupling. It is a trivial 

exercise in second order perturbation theory to show that for a state with momentum q the 

energy changes by: 

6E, = -^Z 



2 V 2MLuJ (?) E i+p -E q + hu (p) 



■Km j p 



*„,2 

In 



AMhce 



2a (m*c s + fiq) 



2a {m*c s — hq) 



(5.6) 



where the integral was cut-off at A w 1/a. For hq m*c s , this simplifies to 



5E In ( h I l) ^ {H + Aam * Cs) hV (5 7) 

9 2Mhc s \2am*c s j 4aMm*c 3 s {% + 2am*c s f 



By comparing (5.7) with (|5.4j ) we conclude that the energy of the carriers can be written as: 



A K h 2 q* 

^-~ + w (5 - 8) 

The Mott gap, Ar, and the effective mass of the carriers, tur, have been renormalized by 
the polaronic effect and are given by: 
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7rm *7p ( h 

Afi - AM "2M^: ln 2^: + 



m R = m 



1 + 



h(h + 4am* c s 



8m*Mc^a 2 (nvp + g 4 ) 2 (h + 2am e c s ) 2 



(5.9) 



and we conclude that the Mott gap is substantially reduced and the mass of the carriers 
increased because of the polaronic effect. 



VI. CONCLUSIONS 

In this paper we have studied the problem of charge transport via non-linear excitation 
in a band and Mott insulators. The insulating behavior in band insulator is due to the 
periodicity of the lattice and the electron-ion interaction opens a gap in the charge spectrum. 
Mott insulators, on the other hand, are dominated by electron-electron interactions and the 
gap is produced by correlation effects when the electronic density is commensurate with the 
lattice. We argue that in band insulators the main propagating excitation is a polaron that 
is described by an electron dressed by the polarization cloud. In the Mott insulator, because 
of the strong electron-electron interactions, the main non-linear excitation is the polaronic 
soliton, a non-linear excitation of the Luttinger liquid that is dressed by a polarization cloud. 

We have shown that in a band insulator the doped electron (or hole) is strongly dressed 
by the internal degrees of freedom in the system that we described in terms of phonons and 
pseudo-spins. We show that the polaron can propagate freely: the lattice and the two-level 
systems follow its motion in the system. Moreover, we have shown that the pseudo-spins 
lead to an extra stabilization, relative to the phonon problem, by reducing its energy and 
size. 

In the case of the Mott insulator we have shown via a bosonization calculation that 
the non-linear excitations can be thought of as a soliton comprising an electron dressed by 
the phonons and pseudospins. We term this new excitation the polaronic soliton. Many 
properties of these excitation are similar to the polaron problem. Namely, the polaronic 
solitons are very stable and their energy is reduced relative to the problem without phonons 
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and pseudo-spins. 

The great qualitative difference between a Mott insulator and a band insulator is in 
the spin spectrum. Spin excitations in a Mott insulator are decoupled from the charge 
excitations (spin-charge separation) and remain gapless even when the Mott gap opens. 
Thus Mott insulators have strong antiferromagnetic correlations. This should be contrasted 
with the case of a band insulator that has both a charge and a spin gap: in order to excite 
an electron from the valence band to the conduction band we have to pay an energy cost 
to unpair two electrons with opposite spins in the valence band. In order to illustrate this 
point let us consider a ID system that is at half-filling with a weak periodic potential, U(x), 
added such that this potential doubles the unit cell. Because the system is half-filled, a 
potential that doubles the unit cell has to open a gap at the Fermi energy. (This mechanism 



is equivalent to the Peierls mechanism due to lattice distortions [|14] , [Lq| ) . Let us consider 
this process of opening a gap using the bosonization technique. We consider the problem 
before the gap opens. We can bosonize the system as explained in Section |. The coupling 
between the weak periodic potential and the electrons can be written as: 

Hu = [ dx^U(x)tlji(x)Mx), (6-1) 

where, because of the periodicity, 

U(x) = U cos{7ix/a) . (6.2) 

We now expand the electron operators as in (|2.3|) in terms of right and left moving electrons 
and find that in Hjj there are terms that oscillate with cos(7nr/a) and cos(7rx/a ± 2kpx). 
Since kp = n /(2a) the terms with wave-vector n/a ± 2kp do not oscillate. All the other 
terms vary very rapidly in space and can be disregarded. In summary, we can rewrite Hjj 



as 



H v = U dx^2i! j Ra {x)ip Lta {x) + h.c. 

— — ■ I dx cos(2y / 7r0 o -(x)) 
na J ~ 

2U t 

= — - / dxcos(V2^$Jx)) cos(v / 2T$ CT (x)) , (6.3) 
7ia J 



23 



where we have used the bosonization rules (|2.4|) and ( |2.5| ). Notice that ( |6.3|) describes the 
backscattering of the electrons by the periodic potential of the lattice. In fact the condition 
that 7f/a = 2kf is nothing but the Bragg condition for this particular problem. We have thus 
learned in the bosonized language a very simple fact about electron scattering in periodic 
potentials. Note, however, that in the bosonic language the extra periodic potential is a 
highly non-linear operator in terms of the bosons. In practice, this term is not easy to treat 
exactly. Let us consider, however, a perturbative RG calculation like the one explained in 
Section |. It is easy to show that: 

dU („ K n K n 



which demonstrates that the periodic potential is a relevant perturbation of the Luttinger 
liquid when K p + K a < 4. Notice that when this happens (as we explained in Section [D 
Hjj becomes relevant and the fields $ p and $ CT are pinned at the minimum of the potential 
energy. In this case, both a charge gap and a spin gap open in the spectrum! The size of 
such gaps, as in fl2.19| ), can be obtained from the RG ( |B.4[ ) and read: 



A„ = ^ (6.5) 

This result shows that a gap in a correlated band insulator depends on the strength of the 
electron-electron interactions. In particular, in the non-interacting limit when K p = K a = 1 
we see that A° M = Uq, as expected from elementary solid state physics. It is clear that in the 
case of a band insulator the spin gap has nothing to do with superconducting fluctuations 
but rather with the fact that the band is full and therefore that the electronic shells are 
closed. 

The RG flow for the periodic potential ( |6.4j ) should be contrasted with the RG flow 
for the Umklapp term, Eq. fl2.17|) , which opens the Mott gap, and is relevant for K p < 1. 



Notice that there is a competition between the Umklapp term and the periodic potential: 
for < K p < K a /3 the Umklapp term diverges faster than the periodic potential while 
for K a /3 < K p < 4 — K a the electron-ion interaction has a stronger divergence. (For 

24 



K p > 4 — K a both operators are irrelevant). Since K p decreases with the strength of the 
electron-electron interaction it becomes obvious that for strong electron-electron interaction 
the Umklapp term is more important than the periodic potential and a Mott gap opens 
in the system. The spin spectrum is gapless while the charge spectrum has a gap. On the 
other hand, if the electron-electron interactions are weaker, the periodic potential dominates 
and a band gap opens in both the spin and the charge spectrum. It is clear that we can 
differentiate the band gap and the Mott gap by examining the spin spectrum. The critical 
value of interactions, U c and V c in (|2.2j ), above which a Mott gap appears in a band insulator 
occurs when K p = K a /3. Using ( p.9|) we find 

U c -h = ^, (6.6) 

which gives the critical line in the U x V plane above which one finds a Mott insulator and 
below which we have a band insulator. 

The applications of the formalist developed here to systems like DNA depends very much 
on the actual values of U and V. If DNA is better described by a band insulator then the 
formalism developed in section [TIT] is more appropriate and the Davydov's equations ( J3.10 ) 



should apply. In this case the elementary excitations should be polarons. However, if the 
electron-electron interactions are strong, as we conjecture, then the polaronic soliton of the 



Mott insulator, as described in section [TV], is the best description. While numerical simu- 
lations can provide useful information on order of magnitude of the electronic parameters, 
only well-controlled experiments in periodic DNA sequences can actually give the final word 
on the insulating nature of DNA. As we discussed previously, we believe that the correct 
way to find out about the insulating nature of DNA is by the study of the spin excitation 
spectrum and not by conductivity measurements that are insensitive to the spin degrees of 
freedom. While there are indications of antiferromagnetic fluctuations in DNA that would be 
consistent with the Mott insulator picture ]7[ the experimental picture is far from complete 
and many more experiments are needed. 

In summary, we have shown that non-linear excitations in band and Mott insulators are 
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very similar and have solitonic character. We show that the band gap is strongly renor- 
malized by electron-electron interactions in ID and that a ID system can be driven from 
a band insulator type behavior to a Mott insulator behavior as a function of interactions. 
The correct way to distinguish between these two types of insulator it is not by studying the 
charge degrees of freedom, that propagate in the form of non-linear waves, but to the spin 
degrees of freedom that remain gapless in the Mott case and are gapped in a band insulator. 
We believe that the physics described here can be applied to strongly correlated systems 
like DNA. 
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